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The magneto-electric properties of resonant tunneling double barrier structures using GaMnAs 
for the quantum well is investigated within a self-consistent Green's function approach and a tight- 
binding electronic structure model. The magnetic state of the well is determined self-consistently by 
the tunneling current which controls the hole spin density and, hence, the degree of exchange splitting 
of the subbands inside the well. Prompted by recent experiments we compare model systems of 
increasing defect concentration (substitutional disorder) regarding their I-V curve, magnetic state, 
and spin polarization. We predict that, near resonance, the ferromagnetic order which may be 
present at zero bias in the GaMnAs well tends to be destroyed. Resonance peaks are found to be 
more sensitive to disorder than ferromagnetic ordering and spin polarization of the steady-state 
current. 

PACS numbers: 85.75.Mm, 73. 23. Ad, 73.63.-b, 72.25.Dc 
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I. INTRODUCTION 



The realization of electric control of ferromagnetism 
in nanostructures is of great interest both for spintronic 
device application and for achieving a better understand- 
ing of the physical mechanisms and dynamics underlying 
the formation of ferromagnetic order in Mn doped semi- 
conductors. Dilute magnetic semiconductors (DMS) are 
made magnetic by doping of ZnS-structured semiconduc- 
tors with transition metal elements, which provide local 
magnetic moments arising from open electronic d or / 
shells [3, 0- A prototype is bulk Gai-^Mn^As in which 
Mn residing on the Ga site (MnQa) donates both a hole 
and a local magnetic moment. Mnc a is an at least mod- 
erately deep acceptor and associated levels lie about 100 
meV above the valence band edge ||. In a recent scan- 
ning tunneling microscopy experiment the radius of the 
Mn acceptor wave function has been determined to be 
about 2 nm Q. Due to an antiferromagnetic exchange 
coupling between the itinerant holes and the local Mn 
d-electrons, an effective ferromagnetic ordering among 
the Mn-ions, known as carrier-mediated ferromagnetism, 
can be established [1-0] • Due to the hole-concentration- 
dependent effective exchange field the spin degeneracy 
of the holes is lifted resulting in a self-consistently spin 
polarized hole gas. The critical temperature T c for the 
occurrence of ferromagnetism in bulk GaMnAs is typi- 
cally below ~ 150 K and depends on Mn concentration 
and sample preparation [H, U, Q . 

Probably due to the degree of (unwanted) defects in 
GaMnAs samples depending on growth conditions, ex- 
perimental evidence has led to somewhat conflicting con- 
clusions about the electronic structure in the vicinity of 
the Fermi level. While all experimental studies on ferro- 
magnetic bulk GaMnAs confirm it to be p-type, there is 
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some debate as to the precise position of the Fermi en- 
ergy . Some experiments can be interpreted by placing 
it into the top of a GaAs-like valence band edge which, 
at most, broadened by disorder Others suggest the 
existence of an isolated impurity band (of high "effective 
mass") which forms at Mn concentrations above ~ 1.5% 
leading to a metal-insulator transition in high-quality 
GaMnAs 0, Hi). 

Irrespective of the detailed electronic structure it is 
clear that, due to disorder, the localization length of va- 
lence band eigenstates will tend to decrease as one moves 
from the top of the valence band towards the energy 
gap Near the band edge a coexistence of local- 

ized and delocalized (Bloch-like) eig enstates may be ex- 
pected, similar to amorphous Si [Tl|, LL2|. ^ Mng a is the 
main defect to provide modification of the valence band 
edge from what it is for GaAs, quantization effects can be 
expected for a layer thickness < 3 nm. Recent tunneling 
spectroscopy of GaMnAs quantum well structures has in- 
dicated such effects [9j. However, the signatures in the 
current-voltage characteristics appear to be rather weak 
and no regions of negative differential conductivity due 
to resonances associated with GaMnAs well layers have 
been observed yet, with the exception of an asymmet- 
ric magnetoresistance resonant tunneling structure [l3| . 
This suggests that a significant concentration of defects 
may be present, as it is the case, e.g., in thin layers of 
amorphous Si, in which similarly weak signatures have 
been found [ill, [HI • Disorder on the other hand, as we 
know from bulk Gai-^Mn^As, goes hand in hand with 
ferromagnetic order and, indeed, ferromagnetic behavior 
has been verified experimentally for thin layers of GaM- 
nAs H03. 

Due to spin-selective hole tunneling in and out of a 
ferromagnetic Mn-doped GaAs quantum well, regardless 
whether occurring sequentially or resonantly, one can ex- 
pect hole population and spin polarization to become de- 
pendent upon the bias applied to a double barrier struc- 
ture containing such a well. This mechanism should allow 
an electric control of ferromagnetic order in the quan- 
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turn well and is evaluated in this paper with particu- 
lar attention to the qualitative effect of disorder in the 
Mn doped well region. An important question addressed 
is whether disorder and unwanted defects can suppress 
spin-selective tunneling. Disorder has a strong influence 
on the electronic properties and they, in turn, strongly in- 
fluence (resonant) tunneling. High-quality ferromagnetic 
layers integrated in semiconductor heterostructures, have 
been envisioned to allow spin-dependent carrier transmis- 
sion [l5[ . Magnetic resonant tunneling structures of high 
structural quality may allow the realization of spin valves, 
spin filtering, and spin switching devices, as proposed in 
several studies [H, Il6l - l22l | . 

In this article we investigate spin-selective hole trans- 
port in GaAs/AlGaAs/GaMnAs/AlGaAs/GaAs hetero- 
structures within the limit of moderately thin samples so 
that an effective independent particle model provides a 
good first approximation. We apply a non-equilibrium 
Green's function formalism based on a tight-binding 
Hamiltonian for the electronic structure, including self- 
consistency regarding the charge density and the ex- 
change splitting of the effective potential, as well as 
charge transfer to the contacts. The carriers' Coulomb 
interaction and the exchange coupling with the magnetic 
ions are described within a mean-field picture. Details of 
our model arc exposed in Sect. [TTJ Since disorder seems 
to play a major role in actual samples we study the ef- 
fect of substitutional disorder on the I-V characteristics, 
the ferromagnetic state of the heterostructure, and spin- 
polarization of the current density. Results and relevance 
to experiment are discussed in Sect. IIII1 Summary and 
conclusions are given in Sect. IIV1 

II. SELFCONSISTENT TRANSPORT MODEL 

The basic features of the semiconductor double-barrier 
structure near the top of the valence band edge are 
mapped onto a two-band tight-binding Hamiltonian for 
the heavy holes (J3 = ±3/2) 

H s = J~]ei t<r \i,<T){i,o-\ 

i,cr 

+ ^2ti,tr^\i,o){i + l,a f \+^c, (1) 

i,crcr' 

where £,- )CT is the spin-dependent (0 =[,]=. ±1) onsite 
energy at lattice site i, t^ aa i denotes the hopping-matrix 
between neighboring lattice sites, and h.c. abbreviates 
the Hermitean conjugate term. Spin conserving hopping 
gives a diagonal matrix ij )Cr ,cr' — t8 aa > with the hopping 
parameter t — ~h 2 /(2m*a 2 ) depending on the effective 
mass m* and the lattice spacing a between to neighboring 
lattice sites. The onsite energy 

ei,a = Ui-e<t>- ~A 4 (2) 

includes the intrinsic hole band profile Ui due to the band 
offset between different materials, the electrostatic po- 



tential (f> with e denoting the elementary charge, and the 
local exchange splitting A<. Near the band-edges this 
model is equivalent to an effective-mass model, however, 
it has the advantage that structural imperfections, as well 
as spin-flip processes, can be readily be modeled by vary- 
ing the onsite and hopping energies. A more realistic de- 
scription can be achieved by introducing a larger set of 
orbital basis functions at each lattice site [15|, I23l425j . 

Within a mean-field approach the exchange coupling 
between holes and magnetic impurities can be described 
by two interrelated effective magnetic fields, respectively, 
originating from a nonvanishing mean spin polarization 
of the ions' d-electrons ( S z ) and from the hole spin den- 
sity ( s z) = (if ~ |26l - [28| . The exchange splitting 
of the hole bands is given by 

A(z) = -Jpdn im p(z)(S z )(z) , (3) 

with z denoting the longitudinal (growth) direction of 
the structure, J p d > is the exchange coupling between 
the impurity spin and the carrier spin density (in case 
of GaMnAs p-like holes couple to the d-likc impurity 
electrons), and fij mp (z) is the impurity density profile 
of magnetically active ions. The magnetic order between 
the impurities is mediated by the holes and the effective 
impurity spin polarization depends on the mean hole spin 
polarization via 

<*,> - -SB S (^g^) , (4) 

where, respectively, ks, T, B$ is the Boltzmann con- 
stant, the lattice temperature, and the Brillouin function 
of order S, here with S — 5/2 for the Mn impurity spin. 
Combation of Eq. ^ and Eq. (Q} gives a self-consistent 
effective Hamiltonian for the holes H c g — — aA(z)/2 with 

A(z) = Jpdn im p(z)bB s < — — J- . (5) 

This shows that a manipulation of the hole spin density 
(s z ) by the applied bias is the key to the control of fer- 
romagnetic order in the heterostructure. 

Within a Hartree mean-field picture space-charge ef- 
fects are taken into account self-consistently by calculat- 
ing the electric potential from the Poisson equation, 

_jL e _l = e [Na{z) _ n(z)] f (6) 

where e and N a , respectively, denote the dielectric con- 
stant and the Mnc a density. The local hole density at 
site |i) is computed as 

-« = £E / ^G<(E;ia,ia) , (7) 

fc||,<7 

with A and k\\, respectively, being the in-plane cross sec- 
tional area of the structure and the in-plane momentum. 
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The non-equilibrium "lesser" Green's function G < is cal- 
culated from the equation of motion 



G< = G R E<G A 



(8) 



where G R and G A — [G fl ] + denotes the retarded and 
advanced Green's function, respectively. The scattering 
function £< = T,f + describes particle inflow of from 
the left (I) and right (r) reservoir and [29| 

£,< =/o(£- «,r) (£&.-££.) , (9) 

where fo(x) = [1 + exp^/fceT)] -1 is the Fermi dis- 
tribution function and fii and /Lt r , respectively, denote 
the quasi-Fermi energies in the contacts. The retarded 
and advanced self-energy terms T, R = T, R + T, R and 
T, A = + , respectively, couple of the simulated sys- 
tem region to the left and right contacts, modeled as 
semi-infinite chains, for which analytic expressions are 
available [ID, The retarded Green's function, finally, 
is given by 



G = [E + ir)-H s - £ 



(10) 



A self-consistent calculation of the spin-dependent ef- 
fective (one-particle) potential, based on the band split- 
ting given by Eq. the Poisson equation Eq. (I6I7[) . and 
the kinetic equations Eqs. ([5]) and (ITU1) . also entails an ad- 
justment of the quasi-Fermi energies relative to the band 
edges in the contacts to ensure asymptotic charge neu- 
trality 31]. These loops must be solved self-consistently 
until convergence to a steady-state solution is reached. 
A small external magnetic field is applied initially to aid 
spontaneous symmetry breaking. For the next bias iter- 
ation, the self-consistent solution from the previous bias 
value is used for an initial guess. In most cases a sta- 
ble self-consistent solution is found. However, due to 
the nonrelativistic nature of the model with orbital mo- 
tion decoupled from the spin degree of freedom, "spin 
up" and "spin down" are degenerate. Moreover, there 
are bias regions where no single stable solution exists re- 
garding charge distribution (intrinsic charge bistability) 
and/or magnitude of the exchange interaction (magnetic 
multi-stability) jaU,!!]. 

Having obtained the self-consistent potential pro- 
file the spin-dependent transmission probability T a i a (E) 
from the left to the right reservoir is calculated from spe- 
cial matrix elements of the retarded Green's function 124] 



T i — T 

- 1 - a tr — a 



v r ,*>\G R {E;r<T'M)\' 
v^G^E-MM)? 



(11) 



with Go denoting the free Green's function of the asymp- 
totic region, and w; )Cr and v TtCr , respectively, are the spin- 
dependent group velocities in the leads. G R (E;ra' ,la) 
is computed most conveniently by adding one layer af- 
ter the other which requires solely 2x2 matrix inversions 
within the present two-band model [30]. 



The steady-state current is calculated within station- 
ary scattering theory (a generalized Tsu-Esaki formula), 



Ja'a 

9(E) 



em*k B T f 00 , 



(2 

, , l + exp[(/i; -E)/k B T] , 
l S l + exp[( Mr -£;)/fc B T] f- - 



The applied bias V = (/i; — jj, r )/e is defined as the dif- 
ference in quasi-Fermi levels of the contacts. In such 
a nonequlibrium situation the hole concentration in the 
GaMnAs layer is not determined not simply by its ferro- 
magnetic state, but also depends on other physical quan- 
tities, such as doping level in the contact regions, struc- 
tural parameters of the heterostructure, and applied bias. 



III. RESULTS AND DISCUSSION 

In this section we employ the model outlined above 
and study the influence of disorder on the resonant tun- 
neling characteristics and the magnetic state of a double- 
barrier structure featuring a GaMnAs quantum well. For 
this simulation we use parameters for GaAs and GaM- 
nAs which are we well established in the literature: m* = 
0.4 m , e r = 12.9, V hax = 300 meV, m = fi r = 70 meV, 
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''imp 



= 1 x 10 20 cm^, J pd = 0.15 



d = 20A, w 

cV nm 3 [33], T — 4.2 K, where mo denotes the free 
electron mass, e r is the relative permittivity, Vbar is the 
bare barrier height of AlGaAs relative to GaAs, d and 
w, respectively, are the barrier and quantum well width. 
The background charge N a arises from an assumed 10% 
of the Mn doping rti mp , since GaMnAs is known to be a 
heavily compensated semiconductor [g, 34]. Doping in the 
contacts, experimentally Be may be used, is such that the 
thermal equilibrium position of the Fermi energy /ij = \i r 
lies near the first heavy-hole resonance, thus, providing 
ferromagnetic order at zero bias. 

Since GaMnAs represents a ternary alloy, disorder 
from Mnca substitution must be accounted for to cap- 
ture its electronic structure. In other ternary alloys, such 
as AlGaAs or CdZnSe, disorder effects lead to features, 
such as direct to indirect gap transitions and band bow- 
ing [Hip. 

In addition, unwanted defects, such as antisites, inter- 
stitials, Mn complexes, vacancies, etc, as well as addi- 
tional defects arising from the heterointerfaces must be 
expected or have been observed @, 3 ■ Clearly, the com- 
plexity of the material is such that a realistic predictive 
electronic structure calculation from first principles is 
practically impossible. Hence we must take a completely 
different approach and model disorder phenomenologi- 
cally, taking doping level and basic properties of MnQ a 
for guidance. Here disorder effects in the GaMnAs layer 
are modeled by performing a configurational average over 
structures with randomly selected onsite and hopping 
matrix elements of the tight-binding Hamiltonian in the 
Mn doped region. We choose a 5% MnQ a concentration 
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FIG. 1: Averaged total density of states near the center of 
the GaMnAs well region for zero bias and maximum hopping 
disorder (a t = 20%). For valence band states (E > 0) integra- 
tion over in-plane momentum is taken into account and the re- 
sult is normalized to the 2D-density of states Do = m* /(nti 2 ). 



in the well and model random (uncorrelated) substitu- 
tional disorder. If a Mn ion is present at a given lattice 
site in the well the onsite energy is shifted according to a 
Gaussian distribution around a mean onsite energy-shift 
of 80 meV and a standard deviation of 20 meV, which 
are reasonable values according to data available for the 
isolated Mnc a acceptor, as well as recent experimental 
results for Gai.^Mn^As for x 2-6% [4, 9j. The hopping 
matrix element is sampled according to a Gaussian with 
5, 10, and 20% standard deviation (at) of its bulk value 
t. This increase is in hopping matrix variation is used to 
simulate an increasing level of defect concentration, for 
fixed MnQa concentration. For each such randomly se- 
lected Hamiltonian the transport problem is solved self- 
consistently and the I-V curve is computed. Final re- 
sults for I-V curve, magnitude of magnetization and spin 
current polarization, etc., are obtained by averaging over 
individual results obtained for these configurations. Typ- 
ically 300 configurations are used to perform this average. 
In this simulation we keep constant the number of active 
Mn spins so that the formation and degree of ferromag- 
netic order at given bias is determined self-consistently 
from the hole spin polarization in the well region. 

Figure [T] shows the averaged density of states near 
the band edge in the Gai-^Mn^As well region for zero 
bias and maximum degree of hopping disorder modeled 
here (a t = 20%). It is seen that the main effect of our 
model for substitutional disorder leads to to a broad- 
ened isolated impurity state above the valence band edge, 
whereby the latter clearly displays the steps characteris- 
tic for 2d quantization effects in spite of disorder (note 
the logarithmic scale!). The defect levels do not provide a 
genuine band since they are confined to the Gai-^Mn^As 
well region. These states may trap holes and contribute 
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FIG. 2: (Color online) Spin-dependent transmission proba- 
bility of the double barrier structure at zero bias with and 
without disorder (a t = 5%). 



to ionized impurity scattering but are not actively in- 
volved in tunneling here since they lie below both quasi 
Fermi levels at up to moderate bias. The transmission 
probability versus energy of the incident holes at zero 
bias is displayed in Fig. [2] In absence of disorder, the 
GaMnAs top valence band structure is modeled as that 
of GaAs plus a self-consistent exchange and the reso- 
nances, indicated as dashed lines in Fig. [2j are spin-split 
by about 30 meV. Taking into account disorder, see solid 
curves for a t = 5% in Fig. [5]) leads to spectral broad- 
ening and a shift of the resonances deeper into the va- 
lence band (anti-bonding effect) of the resonances. An 
increase in overlap of the transmission peaks for spin- 
up and spin-down holes under disorder is particularly 
pronounced for the first heavy-hole resonance since it is 
most sensitive to potential fluctuations. For low contact 
temperatures, however, this effective spin splitting en- 
sures spin-dependent tunneling rates even at a level of 
disorder where the non-monotonic increase of the current 
with applied bias is practically lost, as shown below. 

Contributions from the light-hole band lead to ad- 
ditional resonances. Their inclusion would call for a 
higher-dimensional tight-binding model which properly 
captures nonparabolicity parallel to the heterointerfaces 
and goes beyond the scope of this paper. We just point 
out that, for the present structure a light -hole-band res- 
onance would be expected somewhere between the first 
two heavy-hole-associated resonances, contributing to a 
further masking of negative differential conductance [9] . 

Discontinuous second derivatives of the IV-curve ob- 
tained from tunneling spectroscopy have been attributed 
to quantization effects in the GaMnAs quantum well. Re- 
gions of negative differential conductivity, however, have 
not been observed directly in the IV-curve @ • This sug- 
gests that disorder may play a considerable role, similar 
to amorphous silicon quantum wells for which weak sig- 
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FIG. 3: (Color online) /^-characteristics of a magnetic dou- 
ble barrier structure due to heavy-hole associated bands for 
different degrees of disorder, i.e., standard deviations a t of 
the hopping matrix elements, as explained in the text. The 
inset shows the IV-curve in the voltage range in which the 
transport takes place via the first heavy-hole subband. 
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FIG. 5: Logarithmic local density of states (LDOS) as a func- 
tion of energy at the bias V — 0.04 V (before the first current 
maximum). The self-consistent band profile is indicated by 
the solid line. 
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FIG. 6: Logarithmic local density of states (LDOS) as a func- 
tion of energy at the bias V = 0.052 V (at the first current 
peak). The self-consistent band profile is indicated by the 
solid line. 



FIG. 4: (Color online) The configuration averaged spin split- 
ting |A| in the quantum well as a function of the applied bias 
for different degrees of disorder. 



natures of resonances have been predicted and observed 
[III [r3 | . The current- voltage I-V characteristics for an 
increasing degree of disorder is plotted in Fig. [3] Our 
model reveals that the first region of negative differen- 
tial resistance, corresponding to the first heavy hole res- 
onance, does not disappear until considerable hopping 
disorder of about 10%, modeled as variance in t, is as- 
sumed. This relatively high value needed to flattened the 
IV-curve in our simulations suggests that in real samples 
defects other than MnQ a , such as interstitials, antisite de- 
fects or voids, may play a considerable role for blurring 
the resonances. 



In Fig. 3] the average exchange band spin splitting |A| 
in the quantum well is plotted versus applied bias, show- 
ing that the degree of ferromagnetic order in this struc- 
ture can be controlled by the applied bias. At the po- 
sitions of the current maxima (compare to Fig. [3]) the 
quantum well is populated by the contacts with unpolar- 
ized holes, resulting in a breakdown of the ferromagnetic 
order. The degree of bias control is moderated for in- 
creasing defect concentrations. However, our averaging 
procedure may exaggerate this reduction since it corre- 
sponds to a physical situation where, in plane, there is 
a coexistence of (uncorrelated) ferromagnetic and non- 
ferromagnetic domains. If on the other hand, in-plane 
ferromagnetic order is established or destroyed coherently 
(in correlated fashion), the degree of bias control is un- 
derestimated by our sampling procedure. This makes us 
suggest to perform spin-sensitive tunneling spectroscopy, 
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FIG. 7: Logarithmic local density of states (LDOS) as a func- 
tion of energy at the bias V = 0.06 V (at off-resonance con- 
ditions). The self-consistent band profile is indicated by the 
solid line. 
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FIG. 8: (Color online) Averaged current spin polarization \Pj\ 
versus applied bias V for different degrees of disorder. 



ideally in presence of an external magnetic field, since the 
current spin polarization is predicted to be a more sen- 
sitive signature to determine the degree of ferromagnetic 
order in the sample than resonances in the I-V curve. In 
fact the latter are predicted to occur at a bias when the 
sample is in the nonmagnetic state only. 

To illustrate the changes in the magnetic state of the 
quantum well near the first current peak at V = 0.052 V, 
we plot the local density of states (LDOS) closely be- 
low resonance (V = 0.04 V, Fig. [5]), at resonance (V = 
0.052 V, Fig. EJ), and above resonance (V = 0.06 V, 
Fig.[7]). At zero bias (not shown in Fig. [5]) ferromagnetic 
order is present, placing the Fermi energy between the 
spin split subband edges. The spin up level is found at a 
lower energy than the spin down level. (Note that we use 
an inverted energy scheme for the valence band.) This 
situation is achieved by matching the contact doping level 



to the lowest (heavy-hole) resonance in the heterostruc- 
ture. As bias is increased both spin-up and spin-down 
subband become accessible from the emitter side, how- 
ever, only the lower spin-up subband is accessible from 
the collector side. This maintains hole spin polarization, 
however, it decreases with increasing bias, leading to a 
continuous decrease in A shown in Fig.HJ While hole spin 
polarization decreases, hole density increases in the well 
region partially screening the applied bias. As the ap- 
plied bias is increased resonance is reached. Hole charge 
and exchange splitting in the well are no longer suffi- 
cient to keep the spin-down subband above the collector 
quasi-Fermi level and ferromagnetic order collapses, as 
shown in Fig. |U Both subbands are equally flooded from 
both contacts leading to an unpolarized hole gas in the 
well and, hence, to a destruction of ferromagnetic order. 
Note the change in the effective potential profile in the 
well from convex up in Fig. [5]to almost linear in Fig. |6l If 
the bias is increased further the levels are pushed closely 
below the emitter valence band edge and resonance to 
the emitter is suppressed. However, any small pertur- 
bation triggers the system into a ferromagnetic state, in 
which the spin up level preferrentially is filled from the 
collector side (see Fig. [7]). Thus, inspection of the LDOS- 
plots shows that the ferromagnetic state in the quantum 
well is determined by the relative position of the well 
subband edges and the contact quasi-Fermi levels. An 
appropriate tailoring of these levels allows to change the 
hole gas polarization under bias and thereby provides 
electrical control of the ferromagnetic state in the well. 
In our case we consider only a two-terminal configuration 
with source and drain contacts but the use of additional 
gates in transverse direction (multi-terminal configura- 
tions) provides an additional control knob to move the 
subbands. Although these structures are very difficult 
to realize in practice, they have been studied in a recent 
experiment |22j |. 

The change of the magnetic state in the well is directly 
reflected in the (collector) current spin polarization, as 
shown in Fig. [5J If the entire quantum well becomes 
nonmagnetic, the spin density in the well vanishes and 
the collector current becomes unpolarized. A polarized 
current indicates exchange-split subbands in the well, 
demonstrating ferromagnetism in the well. The exper- 
imental probing of the current spin polarization at the 
collector side therefore would give important additional 
information to confirm the interpretation of recent ex- 
periments regarding size quantization effects in GaMnAs 
quantum wells [9] . 



IV. CONCLUSIONS AND OUTLOOK 

In summary, we have used a steady-state transport 
model to investigate the role of structural disorder on 
the interplay of ferromagnetic order and resonant tunnel- 
ing in double barrier structures with a GaMnAs quan- 
tum well. Ferromagnetic exchange, as well as the hole 
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Coulomb interaction are treated within a self-consistent 
mean field approximation. Disorder effects are modeled 
by random variation of onsite and hopping matrix ele- 
ments of the tight-binding Hamiltonian according to ba- 
sic experimental findings on the MiiQa acceptor. In this 
work we have modeled an electronic structure with an 
isolated impurity band, as supported by most of the re- 
cent experiments. 

For samples which, at zero bias, exhibit ferromagnetic 
order in the GaAsMn well we predict that ferromagnetic 
order is destroyed under bias near (the first heavy-hole) 
resonance. While at resonance the well region is flooded 
by holes these are overall unpolarized thus prohibiting 
the communication of ferromagnetic ordering amongst 
MnQa sites in the well. Within our model we thus are 
able to provide a possible explanation for the absence 
of exchange splitting near resonances, as observed in re- 
cent tunneling spectroscopy measurements on thin GaM- 
nAs layers [9(. Furthermore we find that an experimen- 
tal investigation of the spin polarization of the collec- 
tor current can give information about the presence of 
magnetism in the quantum well. Such a measurement 
is more revealing regarding the ferromagnetic state than 



the search for resonances in the I-V curve, since the for- 
mer is more robust against disorder than the latter. We 
find that disorder tends to suppress negative differential 
resistance region in the IV-curve. Although our model 
is merely qualitative it indicates that substitution disor- 
der from Mn doping alone is not sufficient to explain the 
absence negative differential conductivity in experiment. 

There are a number of open question which should be 
addressed in the future. One pertains to the correlation 
length of ferromagnetic order parallel to the heteroint- 
erface. It has direct influence on the spin- valve action 
and its bias-control of ferromagnetic heterostructures. 
Furthermore, one may ask to what extent magneto- 
transport experiments can distinguish between the two 
main electronic structure models (with and without iso- 
lated impurity band) discussed in the literature. 
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